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Abstract. We study the first microlensing event toward 
the Carina spiral arm discovered by the OGLE collabo- 
ration. We demonstrate that this long duration event ex- 
hibits strong parallax signatures. Additional information 
from the parallax effect allows us to determine the lens 
transverse velocity projected onto the Sun-source line to 
be ^ 145 km s~^. We also estimate the optical depth, event 
rate and duration distribution for microlensing toward Ca- 
rina. We show that this event is broadly consistent with 
these predictions. 
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1. Introduction 

Gravitational microlensing was originally proposed as a 
method of detecting compact dark matter objects in the 
Galactic halo (Paczyhski 1986). However, it also turned 
out to be a powerful method to study Galactic struc- 
ture, mass functions of stars and extrasolar planetary sys- 
tems (for a review, see Paczyhski 1996). Earlier microlens- 
ing targets include the Galactic bulge, LMC, SMC and 
M31. Recently, the EROS and OGLE collaborations have 
started to monitor spiral arms. The EROS collaboration 
(Derue et al. 1999) has announced the discovery of three 
microlensing events toward two spiral arms. In this pa- 
per, we study the first spiral arm microlensing event dis- 
covered by the OGLE II experiment (Udalski, Kubiak & 
Szymahski 1997). This ongoing event, OGLE-1999-CAR- 
1, was discovered in real-time toward the Carina arm by 
the OGLE early-warning system (Udalski et al. 1994). We 
show that this is a unique event that exhibits strong par- 
allax effects. Such events were predicted by Refsdal (1966) 
and Gould (1992). The first case was reported by the MA- 
CHO collaboration toward the Galactic bulge (Alcock et 
al. 1995). OGLE-1999-CAR-1 is the first parallax event 
discovered toward any spiral arm. The outline of the pa- 
per is as follows. In Sect. 2, we briefly describe the ob- 



servational data that we use. In Sect. 3, we present two 
different fits for the light curve, with and without consid- 
ering the effect of parallax and blending. In Sect. 4, we 
calculate the expected optical depth, event rate and du- 
ration distribution toward Carina. Finally in Sect. 5, we 
summarize and discuss the implications of our results. 

2. Observational Data 

The observational data were collected by the OGLE col- 
laboration in their second phase of microlensing search 
(Udalski et al. 1997). The search was done with the 1.3- 
m Warsaw telescope at the Las Campanas Observatory, 
Chile which is operated by the Carnegie Institution of 
Washington. The targets of the OGLE II experiment in- 
clude the LMC, SMC, Galactic bulge and spiral arms. 
All the collected data were reduced and calibrated to the 
standard system. For more details on the instrumentation 
setup and data reduction, see Udalski et al. (1997, 1998). 

The event, OGLE-1999-CAR-1, was detected and an- 
nounced in real-time on Feb. 19, 1999 by the OGLE col- 
laboration!^. Its equatorial coordinates are a=ll:07:26.72, 
(5=-61:22:30.6 (J2000), which corresponds to Galactic co- 
ordinates I = 290°.8,6 = -0°.98. The ecliptic coordi- 
nates of the lensed star are A = 331°. 9, /3 = — 58°.l (e.g., 
Lang 1980). Finding chart and Lband dataof OGLE-1999- 
CAR-1 are available at the above web site and the V-band 
data were made available to us by Dr. Andrzej Udalski. 
In total, there are 416 I-band and 85 V-band data points. 
The baseline magnitudes of the object 0, in the standard 
Johnson V and Cousins I-band, are V — 19.66, / ~ 18.01. 

3. Models 

In this section, we will fit the OGLE V and I-band light 
curves simultaneously with theoretical models. We start 
with the simple standard model, and then consider a fit 
that takes into account both parallax and blending. 
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Most microlensing light curves are well described by 
the standard form (e.g., Paczyhski 1986): 



L^U^ +4' 



u{t) = ^ul + w-^{t), 



(1) 



where uq is the impact parameter (in units of the Einstein 
radius) and 



w(t) = ^ tE = RE/vt 

tE 



(2) 



with io being the time of the closest approach (maximum 
magnification), Re the Einstein radius, vt the lens trans- 
verse velocity relative to the observer-source line of sight, 
and Ie the Einstein radius crossing time. The Einstein 
radius is defined as 



Re 



AGMDsx{l - x) 



(3) 



where M is the lens mass, Ds the distance to the source 
and X = D^/Ds is the ratio of the distance to the lens and 
the distance to the source. 

To fit both the I-band and V-band data with the stan- 
dard model, we need five parameters, namely. 



uo,to,tE,mifi,my^o. 



(4) 



Best-fit parameters are found by minimizing the usual 
using the MINUIT program in the CERN library and are 
tabulated in Table 1. The resulting is 893.9 for 496 de- 
grees of freedom. For convenience, we divide the data into 
one 'unlensed' part and one 'lensed' part; the former has 
t = J.D.- 2450000 < 1150d and the latter has t > 1150d. 
For the standard fit, the lensed part has = 407.9 for 
161 data points, and the unlensed part has — 486.0 for 
340 data points; the somewhat high for this part may 
be due to contaminations of nearby bright stars (as can 
be seen in the finding chart), particularly at poor seeing 
conditions. The per degree of freedom for the 'lensed' 
part is about 2.5, indicating that the fit is not satisfac- 
tory. This can also be seen from Fig. 1, where we show 
the model light curve together with the data points. As 
can be seen, the observed values are consistently brighter 
than the predicted ones for t > 1325 d in the I-band. Fur- 
ther, the prediction is fainter by about 0.05 magnitude at 
the peak in the V-band. We show next that both incon- 
sistencies can be removed by incorporating parallax effect 
and blending. 

To take into account the Earth motion around the Sun, 
we have to modify the expression for u{t) in eq. (|l|). This 
modification, to the first order of the Earth's orbital ec- 
centricity (e = 0.017), is given by Alcock et al. (1995) and 
Dominik (1998): 

u^(t) = ul + w{t)'^ +7^1^ sin^tp 

+ 2f® sin -0 [w{t) sin 9 + uq cos 6] + sin^ (3 cos^ -0 
-|- 2r Q sin (3 cos tp [w{t) cos 9 — Uq sin 6] , (5) 



where 9 is the angle between Vt and the line formed by 
the north ecliptic axis projected onto the lens plane, uq 
is now more appropriately the minimum distance between 
the lens and the Sun-source line. The expression of and 
are given by 



and 



lAU 

V tE 



{l-ecospoit-tp)]}, 



ip = -(/) + flait - tp) + 2esin[nQ{t ~ tp)], 



(6) 



(7) 



where tp is the time of perihelion, v — Vt/{1 ~ x) is the 
transverse speed of the lens projected to the solar position, 
f2o — 27r/yr, and 4> is the longitude measured in the eclip- 
tic plane from the perihelion toward the Earth's motion; 
this is given in the appendix of Dominik (1998), 



X + n + (p, 



71 



(8) 



where cpj is the longitude of the vernal equinox measured 
from the perihelion, (p^ — 1.33 (rad), and the Julian day 
for Perihelion is tp = 2451181.57; the readers are referred 
to the The Astronomical Almanac (1999) for the relevant 
data. Note that the inclusion of the parallax effect intro- 
duces two more parameters, v and 9. 

The two-color light curves show that the lensed object 
became bluer by ~ 0.05 mag at the peak of magnification; 
such chromaticity is easily produced by blending. The ad- 
ditional source of light may be from the lens itself and/or 
it can come from another star which lies in the seeing disk 
of the lensed star by chance alignment. When blending is 
present, the observed magnification is given by 



A = h + f^)A{t), l^I,V. 



(9) 



To model the blending in two colors, we need two fur- 
ther parameters - the fraction of light contributed by the 
unlensed component in I and V, // and /y, at the base- 
line. Therefore, a fit that takes into account both parallax 
and blending effects requires 9 parameters: uq, ^Oi ^e, "^1,0, 
mv,o, v,0,fi, and fv- 

The best-fit parameters for this model are given in Ta- 
ble 1. Compared with the standard fit, the x^ is reduced 
from 893.9 to 640.8. The reduction in the lensed part is 
dramatic: the drops from 407.9 to 177.7 for 161 data 
points. The x^ for the unlensed part is 463.2 (as compared 
to 486.0 for the standard fit) for 340 data points. The x^ 
per degree of freedom is satisfactory. The predicted light 
curve (solid line in Fig. 1) matches the observed data both 
in the I-band and V-band. From Table 1, the blending 
fractions in the V and I bands are not well constrained, 
// 0.36±0.2 and fv = 0.29±0.2. The differential blend- 
ing, however, is reasonably constrained, fv/fi = 0-^^-o.27 
due to the observed differential magnification (0.05 mag) 
between the I-band and V-band. The projected lens ve- 
locity is well constrained while its direction has somewhat 
larger errors. For completeness, wc mention that the best 



Mao: An Ongoing OGLE Parallax Microlensing Event Toward Carina 



3 



400 600 800 1000 400 600 800 1000 




1200 1400 1200 1400 

J. D. -2450000 J. D. -2450000 

Fig. 1. The Tband (left) and V-band (right) hght curves observed by the OGLE coUaboration are shown. In each 
panel, the constant part of hght curve is shown at the top with their amphtudes shifted by 1.5 magnitude and with 
their time intervals (from 400 d to 1100 d) labelled at the top axis. The dotted line indicates the best-fit standard 
model (eq. [^) , while the solid line is for the best-fit model that takes into account both parallax (eq. ||) and blending 
(eq. Fit parameters are given in Table 1. 

Table 1. The best standard model (first row) and the best parallax model with blending (second row) for OGLE- 
1999-CAR-l. 



to tE y-o mifi rnvfi 6 v fi fv 

S 1284.2l:«;i 89.8lg;g 0.35lg:°i 18.0lt°;|;^ 19.62t°°^ — — — — 

P 1288.8tl:^ nS.ltZ^ 0.24t°o°l l8-Olt°of2 19.66t°;EJ^ -O-Q'^^ofs WS.Stig;" 0.36t°J 0.29t»i 640.9 



model that accounts for blending but not parallax has 
= 791.8 while the best model that accounts for parallax 
but not blending has = 682.2. Hence the parallax effect 
reduces the much more effectively than blending. This 
can be easily understood since the observed light curve is 
asymmetric, which cannot be produced by blending. 



4. Optical Depth and Event Rate Toward Carina 

To put OGLE-1999-CAR-1 into context, in this section we 
estimate the optical depth, event rate and event duration 
distribution toward Carina. These can be compared with 
future observations when more events become available. 

One major uncertainty for microlensing toward spiral 
arms such as Carina is that we do not know the distance 
to the sources. Several molecular clouds at the same di- 
rection have distances between 6.8 kpc and 8.7 kpc (see 
Table 1 and Fig. 4. in Grabelsky et al. 1988). We adopt 
Ds = 6.8 kpc as the distance to the sources. This as- 
sumption is also consistent with the lensed star being 
a main-sequence star as required by its position in the 
color-magnitude diagram (Udalski 1999, private commu- 
nication); we return to this issue in the discussion. 



Since the direction toward the lensed star is nearly in 
the Galactic plane (6 = — 0°.98) and far away from the 
Galactic center, we assume that the lenses are entirely 
contributed by disk stars, and their density profile follows 
the standard double exponential disk distribution 

p(i?,.) = ^exp(-^)exp(-M), (10) 

where E — SOA/q/pc^, r is the lens distance to the Galac- 
tic center, tq = 8.5 kpc, disk scale-length r^i = 3.5 kpc and 
scale-height zq — 0.325 kpc. The optical depth can then 
be obtained (e.g., eq. 9 in Paczynski 1986): 

r = 3.4xl0"^ (11) 

independent of the lens mass function and kinematics. 

To estimate the event rate and event duration distri- 
bution, we have to make assumptions about the lens kine- 
matics and mass function. The motion of lenses can be di- 
vided into an overall Galactic rotation of 220 km s^^ and a 
random motion. We assume that the random component 
follows Gaussian distributions with velocity dispersions 
of an — 4Oknis~^,(T0 = 30kms^^,(Tz — 20kms~^ (cL 
Derue et al. 1999). The motion of the Sun relative to the 
Local Standard of Rest is taken as vr^q = 9 kms~^, "vg^Q = 
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Fig. 2. The left panel shows the predicted event rate 
distribution (in units of per million stars per year) 
as a function of duration toward Carina. The short 
dashed line is for a (5-mass function of IMq, the long 
dashed line is for a Salpeter mass function, n(M)dM cx 
M~'^-^^dM and the solid line is for a disk mass func- 
tion, n{M)dM oc M^^-^^dM, determined from HST star 
counts. The shaded region indicates the la range in time- 
scale for the best parallax fit with blending. The right 
panel shows the likelihood function and lens mass as a 
function of the lens distance, respectively. 

llkms~^,Vz^Q — 16kms^^. We study three mass func- 
tions: 

1. n{M)dM cx 5{M - lMQ)dM, 

2. n{M)dM cx M-^-^^dM, O.OSM© <M < IMq, 

3. n{M)dM cx M-^-^^dM, O.IM© <M < 0.6Mq. 
Note that the second is a Salpeter mass function while 

the third describes the local disk mass function determined 
using HST star counts (Gould ct al. 1997). The predicted 
event duration distributions for these three mass functions 
are shown in the left panel of Fig. 2. It is clear that there 
is a tail toward long durations for all three distributions. 
The probabilities of having longer than 100 d are re- 
spectively, ~ 50%, 10% and 20%; so the observed long 
duration is not statistically rare. The total event rates per 
million stars per year are found to be 



r = 0.64, 1.4, 1.1, 



(12) 



for the three mass functions, respectively. The predicted 
duration distribution and event rate are sensitive to the 
assumed velocity dispersions. For example, if we adopt 
<^R = <^e — <^z — 20kms~^ (as in Kiraga & Paczyhski 
1994), then the events will be on average 30% longer and 
the event rate decreases by about 30%. 

5. Discussion 

We have shown that the microlensing event, OGLE-1999- 
CAR-1, has a light curve shape that is significantly modi- 
fied by the earth motion around the Sun. This event is still 



ongoing at the time of writing (Aug. 12, 1999); later evolu- 
tions of the light curve will test our predictions and reduce 
the uncertainties in the parameters. For comparison, the 
three microlensing events seen by the EROS collaboration 
(Derue et al. 1999) are toward different spiral arms which 
are somewhat closer to the Galactic center direction. The 
three events have ts between 70 d to 100 d, and none of 
the events show parallax effects. 

We have assumed a source distance of 6.8 kpc in the 
previous section. We show now that this is a reasonable 
assumption. From V = 19.66 and / = 18.01, and taking 
into account the blending, we find that the lensed star has 
My = 5.9, Ml = 4.3. Assuming an extinction of Ay = 1.5 
mag (see Fig. 4 in Wramdemark 1980), and Aj/Ay = 
0.482, we obtain the intrinsic magnitude and color My = 
4.4, {V — I)o = 0.75. The star is consistent with being 
a main sequence star (with mass M ^ 1.05Mq), as can 
be seen from the color-magnitude diagram (Udalski 1999, 
private communication) in the Carina region. 

We can combine the expression for v and ts to obtain 
the lens mass as a function of its distance 



M = 



1 ~ X V t 



1.8M 



1 



O" 



(13) 



Using eq. (6) in Alcock et al. (1995), we have calculated 
the likelihood of obtaining the observed transverse veloc- 
ity and direction as a function of the lens distance. The 
result is shown in the right panel of Fig. 2 for the Gould 
et al. (1997) mass function. The lens distance is approxi- 
mately X = 0.5 ± 0.2; we caution that this likelihood func- 
tion is sensitive to the assumed kinematics. For example, 
if we take uji = ag — — 20kms~^, then the best- 
fit lens distance changes to a; = 0.7. It is possible that a 
lens is a (massive) dark lens such as a white dwarf, neu- 
tron star and black hole; in this case the blended hght is 
contributed by another unrelated star. Another possibil- 
ity is that the blending is contributed by the lens itself. 
If the lens is a main sequence star and x « 0.7, which 
implies M « O.SMq, then the lens will contribute about 
the right amount of light to explain the light curves (in- 
cluding the color change in the V and I bands). Further 
high-resolution imaging of the lensed object should allow 
us to differentiate between these two possibilities. 

Acknowledgements. We greatly appreciate the OGLE collab- 
oration for permission to use their data for this publication. 
I also thank Matthias Bartelmann, Martin Dominik, Bohdan 
Paczyhski, Andrzej Udalski, Hans Witt and particularly the 
referee Eric Aubourg for valuable discussions and comments 
on the paper. 

References 

Alcock C, AUsman R.A., Alves D. et al., 1995, ApJ 454 L125 
Derue P., Afonso C, Alard C. et al., 1999, preprint (astro- 



ph/9903209) 



Dominik M., 1998, A&A 329, 361 



Mao: An Ongoing OGLE Parallax Microlensing Event Toward Carina 



5 



Grabolsky D. A., Cohen R. S., Bronfman L., Thaddeus P., 

1988, ApJ 331, 181 
Gould A., 1992, ApJ 392, 442 

Gould A., Bahcall J. N., Flynn C, 1997, ApJ 482, 913 

Kiraga M., Paczynski B., 1994, ApJ 430, LlOl 

Lang K. R., 1980, Astrophysical Formulae (Springer- Verlag: 

Berlin), 504 
Paczynski B., 1986, ApJ 304, 1 
Paczynski B., 1996, ARA&A 34, 419 
Refsdal S., 1966, MNRAS 134, 315 

Udalski A., Kubiak M., Szymanski M., 1997, Acta Astron. 47, 
319 

Udalski A., Szymanski M., Kaluzny J. et al., 1994, Acta As- 
tron. 44, 227. 

Udalski A., Szyrnariski M., Kubiak M. et al., 1998, Acta As- 
tron. 48, 147 
Wramdemark S., 1976, A&AS 23, 231 



